Few-body bound state stability of dipolar molecules in two dimensions 
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Bound structures among dipolar molecules in multilayers are a topic of great interest in the light 
of recent experiments that have demonstrated the feasibility of the setup. While it is known that 
two molecules in two adjacent layers will always bind, larger complexes have only been scarcely 
addressed thus far. Here we prove rigorously that three- and four-body states will never be bound 
when the dipoles are oriented perpendicular to the layers. The technique employed is general and 
can be used for more molecules/layers and other geometries. Our analytical findings are supported 
by numerical calculations for both fermionic and bosonic molecules. Furthermore, we calculate the 
reduction in intralayer repulsion necessary to bind large complexes and estimate the influence of 
bound complexes in systems with many layers. 

PACS numbers: 03.65. Ge, 36.20.-r, 67.85.-d 



I. INTRODUCTION. 

The production and manipulation of ultracold molecu- 
lar gases is a truly riveting physical adventure at the mo- 
ment [I]. Interacting dipolar molecules confined to two 
dimensions are currently attracting substantial attention 
by way of its recent experimental realization and study 
of geometrical effects on chemical reactions [2 . Theo- 
retically, this reduced dimensional setup has spurred im- 
mense interest and many promising proposals for novel 
many-body states exists [Ji, . Many of these studies ignore 
or make rough guesses as to which few-body states are 
allowed in a setup with adjacent two-dimensional (2D) 
planes. An understanding of the basic building blocks is 
obviously crucial when trying to understand the com- 
plicated correlations in the many-body systems. The 
low-energy few-body spectrum with short-range interac- 
tions is well-known |3] and gives rise to the famous Efi- 
mov spectrum in three dimensions (3D) and universal 
trimer states in 2D [51 [7] • However, the dipolar interac- 
tions is anisotropic and long-ranged and it is not clear 
that similar structures emerge. 

A prominent example of the influence of bound states 
on the many-body physics is the BCS-BEC crossover in a 
two-component Fermi gas with attractive short-range in- 
teractions [5]. In the 2D case, the presence of a two-body 
bound dimer state plays a crucial role in the crossover [9] . 
This remains true for dipolar molecules in a bilayer with 
perpendicular polarization [TD], and an effective dimer- 
ization of the system along these lines have been sug- 
gested in the multilayer setup [H]. In the bilayer case, 
this can be infered from the dipolar potential which al- 
ways supports a bound two-body state [12,- 15^. However, 
there is a gap in our current knowledge of the stability 
of few-body states containing more than two molecules. 
In the present letter we give a complete characterization 
of the bound state structures for up to five layers with 
perpendicularly polarized dipoles (see Fig. [T]). 

The most stable three-body trimer system is a chain 
with one molecule in each layer where the pairwise in- 
teractions are attractive at short distances P51. Two 
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FIG. 1: (color online) Dipolar molecules in a multi-layer struc- 
ture with dipoles oriented perpendicular to the planes, a il- 
lustrates the different chains that are the most stable con- 
figurations. Three-body states are shown in 6i and &2, and 
four-body states in Ci, C2, and C3. di shows a non-chain five 
layer configuration, while ^2 is a picture of a weakly bound 
dimer with twice the layer separation. 



molecules in the same layer will repel each other. The 
most interesting trimer is therefore three molecules in 
two layers where two pairs can form bound subsystems 
whereas one pair is entirely repulsive (61 in Fig. [ij . The 
question is whether this trimer is bound with an energy 
below the dimer and a free molecule (62 in Fig- [l}- This 
is a special case in 2D of a repelling pair where each is 
attracted by a third particle. The four-body tetramer 
system with a fourth molecule in the layer next to the 
repelling pair (ci in Fig. [T]) is much closer to stability 
with five bound two-body subsystems and only one re- 
pelling pair. The question here is if that structure has a 
lower energy than a bound string of three molecules and 
a free molecule (c2 in Fig. [T]) where only three bound 
subsystems remain. This of course implies that adding 
more molecules in favorable configurations could produce 
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stable structures. 

In the present letter we first rigourously prove that 
trimer and tetramer states are unstable for identical 
bosons or fermions in the layers. The mathematical 
method is very general and reduces the problem to a 
matter of geometrically identifying the threshold in the 
potential landscape of the systems. It can be applied 
to other systems where repulsive and attractive terms 
compete, like the study of excitons in carbon nanotubes 
[T71 [T5] or in organic semiconductors [121 [20]. A suit- 
able reduction of the in-plane dipolar repulsion leads to 
bound complexes as all interlayer interactions support 
bound states. We expect that this can be achieved in 
different ways, either by applying an electric field gra- 
dient perpendicular to the planes, by having different 
molecules in the system, or by direct manipulation of 
the dipole interaction in each layer through microwave 
transitions [13 [55] . Here we determine the critical repul- 
sive strengths where stable structures emerge analytically 
and support our findings by detailed numerical calcula- 
tions using the stochastic variational method [23|. 



II. FORMULATION. 

The Schrodinger equation for N particles with masses 
m is 

^ t2 ^ \ 



with coordinates fi and energy E^. The interaction be- 
tween particles i and j is XijVij, i.e. 



where d is the interlayer distance and Di is the dipole 
moment of molecule i. For n = we get the intralayer 
two-body l/r^ repulsion, while n > produces the inter- 
layer potentials when the separation is nd. The strengths 
are defined by the product DiDj, and given through the 
dimensionless quantities A^. The Vij{r,n = 0) is mod- 
ified at small distance by the transverse motion. This 
leads to the expression [24 



where U is the confluent hypergeometric Kummer func- 
tion and I is the width of the layer. We have checked that 
our results are insensitive to I in the interval 0.1 < l/d < 
0.2 and the detailed properties at small distance are not 
important. Note that we assume that there is no tunnel- 
ing between the layers, which implies that symmetriza- 
tion is only relevant for particles in the same plane. 



A. Proof of instability. 

Consider first the three-body system with molecules 1 
and 2 in the same plane and molecule 3 in an adjacent 
plane {bi in Fig. [I]). We assume that Di = D2 which 
may be different from D^. We denote the attractive and 
repulsive strengths by Aq = A13 = A23 and Xr = A12. 
Notice that both Aa and A^ are positive. The reason 
we talk of an attractive interaction in this case is that 
the potential always supports a bound state [TT. We 
shall now prove that the system is unbound for A^ — Xa- 
Clearly, the system is bound for a sufficiently small ratio 
Xr/Xa- We can then, for any given Aa, define a critical 
repulsive strength, AJ;''[Aa] (which is thus a function of 
Aq as indicated by the brackets), as the value where the 
energy precisely equals the two-body energy, i?2[Aa], of 
the bound 13 and 23 (interlayer) subsystems, i.e. 

,c.r, 1 _ E2[Xa]-{^\T+XaV^,+XaV23\^) , 

where T is the kinetic energy operator and ^' is the nor- 
malized wave function at the threshold for binding. We 
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We now make re- 



peated use of the general fact that, for a Hamiltonian h 
with eigenstate \4>) and eigenvalue e, i.e. h\4>) = t\4>), the 
partial derivative with respect to paramters that enter h 
(and therefore also enter e) is e' = {(t)\h'\4)). Applying 
this fact for the parameter A^, we obtain 



aA-[Aa] _ ^^-(*|V^13 + V^23|*) 



dXa 



(*|Vi2|*) 



(5) 



Defining the two-body bound state wave function, $, and 
the related kinetic energy operator, T2, we get 



^aPcr 



\a^^^-E2[Xa\ + {^T\^) 



(^ITI*) 



-($|T2|$) 



(6) 



The structure of is the product, ^ = ifj^, where ip 
only depends on the relative coordinate, y, between the 
dimer, $, and the third molecule. This follows from 
the tail of the dipolar potential which is of shorter 
range than r~^, which prohibits long-range correlations 
[1] . The universal limit with a two-body cluster structure 
is therefore approached [55]. We immediately conclude 
that {^\T\^) - ($|r2|$) = {^\Ty\^) is positive, where 
Ty = T —T2 is the kinetic energy operator depending only 
on y. Thus, Fcr > and the relative critical strength in- 
creases with Xa- 

We now compute AJ;'"[Aa]/AQ for Aq, = 00, which is 
equivalent to infinite masses and negligibly small ki- 
netic energies. The critical repulsive strength is found 
when the total potential in all space becomes larger 
than the two-body energy. Thus, the potential surface. 
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Vi3+V23+Vi2Xr^ [Xa]/Xa, depending on two relative coor- 
dinates, should have a minimum less than the minimum 
of the two-body potential in Eq. This geometrical 
condition yields X'^'^[Xa]/\a < 0.373. Since Fcr > 0, this 
is a strict inequality for all finite values of Xa- In partic- 
ular, the system is unbound for all Xa when = A^. 

The most stable four-body system is a chain with one 
molecule per layer (subset of a in Fig. [I]), followed by the 
ci configuration in Fig. [l] Four molecules in two layers 
are much less stable due to the additional repulsion. To 
study the stability when changing the repulsion, wc con- 
sider two molecules in a single layer with D2 — -D3 and 
one molecules in each of the layers above and below it 
with Di = D4, keeping the product -D1-D3 fixed. The 
strengths are A^ = A13 = A12 = A24 = A34, A^ = A23, and 
Ai4 = A^/Ar. The threshold structure, 4*, approached 
when Xr X'^'^[Xa], is three molecules in three layers, 
and a free fourth molecule (c2 in Fig. [T]) . The condition 
of equal energies analogous to Eq. Q becomes 



(«'|T + A-[A,]V^23 



Vii + Aa(Vi3 + ^12 + ^24 + V. 
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= ($|T3 + XaVi3 + A„l^24 + X^I^mI*). 

Similar to the three-body case, the threshold structure is 
\1/ = ■(/)<i), where tp is the wave function depending on the 
relative coordinate, y, between the trimer and the fourth 
molecule. This means that (\l/|Vi4|*) = ($|Vi4|<J>), im- 
plying that the critical value easily is found from Eq. ([t]). 
A completely analogous procedure as above yields 



(8) 



where T3 is the three-body kinetic energy operator. This 
derivative is postive and an upper bound for stability is 
obtained from A^ = 00. The geometries of the potentials 
are such that A^''[Aa]/Aa < 0.75, and all these four-body 
systems in two or three layers are unstable. 

Our results strongly indicate that additional layers 
could lead to stable configurations with more than one 
molecule in each layer. The proofs above hold for bosons 
and distinguishable particles. However, the threshold for 
fermions is always higher and produces less stable sys- 
tems. 



III. NUMERICAL RESULTS. 

We illustrate our findings numerically using the 
stochastic variational method with correlated gaussians 
P3j . Fig. [2] shows two- and three-body energies as func- 
tion of XjfXa with fixed Xa — 10 for the h\ configura- 
tion in Fig. [1] The trimer energies increase with increas- 
ing repulsion, and become unstable at the dimer thresh- 
old. Resonances could exist at higher energies but the 
effective barriers are tiny for these potentials. A large 
width would immediately be accumulated and the reso- 
nance disappears by merging into the two-body contin- 
uum background. The fermion energies exceed the boson 
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2: (color online) Numerical trimer enegies as function of 
t^T/Aa = of the configuration 61 in Fig. [l] calculated 

witii> Aa = 10 for bosons (solid green) and fermions (dashed 
blue^. A horizontal dashed marks the two-body threshold. 
The inset compares the two-layer dimer (solid red) with the 
three-layer trimer as function of Aa (dot-dashed blue). The 
(green) dashed curve is twice the dimer energy demonstrating 
that the three-layer trimer is more stable than two free dimers. 



energies due to the antisymmetry which itself costs en- 
ergy compared to symmetric systems. This is unavoid- 
able in spite of the repulsion because the attraction at 
short distance must be larger to provide a bound state, 
implying that the two molecules in the same layer cannot 
be far apart. The distance from stability is reflected in 
the critical repulsion below which the system is bound. 
The lower curves in Fig. |3] show the numerical results 
for trimers. At large Xa the analytical results 0.373 is 
approached very slowly, while at small Xa only a small 
repulsive perturbation is needed to break the exponen- 
tially weakly bound system. 



Consider now the four-body system in configuration 
ci of Fig. [1] As demonstrated in the inset of Fig. [2j a 
chain of three molecules in three layers is more bound 
that two dimers. This can be understood since it gains 
energy from the attractive interaction of configuration d2 
in Fig.[l] The threshold for binding ci is therefore a chain 
of three molecules and a free molecule, C2, rather than C3. 
The energies for the ci system are shown in Fig. |4] where 
we again see that the threshold is reached for Xr/Xa < 1. 
The relative critical strengths are shown in Fig. [3j where 
the asymptotic strong interaction limit of 0.75 obtained 
from geometry is approached. 
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FIG. 3: (color online) The critical strengths, X^^ /Xa, as func- 
tions of Xa, for configurations 6i (lower curves) and ci (upper 
curves) of Fig. [l] with fermionic (solid) and bosonic (dashed) 
species. The asymptotic values, 0.373 and 0.75, for large Xa 
are slowly approached. 




FIG. 4: (color online) Tetramer energies as function of Xr/Xa 
in configuration ci of Fig. [l] with Xa = 10 (upper) and Xa = 20 
for fermions (dashed purple) and bosons (dot-dashed blue). 
The solid red curve shows the threshold energy (c2 of Fig.[T|. 



A. Structure. 

The structure of trimer and tetramer is contained in 
relative distances within the bound state. In the strong 
dipole limit, all molecules remain close until the repulsion 
is extremely close to the threshold where one molecule 
eventually leaks out through the repulsive barrier. The 
universal limit of (r^)|_B| — > ^ is approached from be- 
low but only for very small energies [111 [25] . For weak 
dipoles, the wave functions are already distributed over 
large regions of space. Our numerical calculations con- 
firm these expectations that are used in the proof of in- 



stability. 

IV. CONCLUSIONS AND PERSPECTIVES. 

We have demonstrated that in two or three layer sys- 
tem with perpendicularly polarized dipolar molecules, no 
three- or four-body bounds states are stable. This was 
shown rigorously through a novel analytical approach 
that predicts geometrically how much the repulsive inter- 
actions in the system must be reduced to stabilize such 
states. We also presented numerical results in support 
of our conclusions and predicting exact critical values 
for when bound state with more than two molecules are 
stable in the setup. Fermionic systems are found to be 
less stable than bosonic for all systems studied. Here we 
have studied the case of perpendicular polarization, but 
we note that similar arguments can be applied for moder- 
ate tilting angles where the in-plane interaction remains 
repulsive on average. 

While we find that chains are the most stable con- 
figurations (a in Fig. [T]), our results show that four-body 
states are much closer to stability than three-body states. 
This implies that an increase in the number of layers can 
lead to bound complexes with more than one molecule 
per layer (as for instance shown in configuration di of 
Fig. [T]). Disregarding these few-body states, both den- 
sity wave and crystalline [57] phases have been pre- 
dicted in the multilayer setup. With a richer spectrum of 
few-body complexes, we expect that more complicated, 
stable many-body configurations are possible when the 
complexes discussed here are used as building blocks. 
Detecting the structures experimentally should be pos- 
sible by several techniques such as RF spectroscopy [28] . 
lattice shaking [29] , or in situ optical detection [SO] |3T] . 

As demonstrated here, one can tune bound state sta- 
bility by reducing the intra-layer repulsion. If an elec- 
tric field gradient is applied that is steep enough to 
make the interaction in adjacent layer different the ef- 
fective repulsion could be reduced and stabilize few- 
body bound states. Alternatively, one can imagine driv- 
ing microwave transitions among rotational states in the 
molecules [21] [22] into states of different effective dipole 
moment for a given external dc field. Through this tech- 
nique it might be possible to even reverse the direction 
of the effective dipole moment, yielding a potential that 
is repulsive at short distance and attractive in the tail. 
This will produce entirely new structures and is an ex- 
citing direction for future studies. 



Acknowledgments 

We thank J. Arlt, J. F. Sherson, and B. Wunsch for 
enlightening discussions. 



[1] K. Winkler et al, Phys. Rev. Lett. 98, 043201 (2007); S. et ai, Science 321, 1062 (2008); K. K. Ni et ai, Science 

Ospelkaus et ai, Nature Phys. 4, 622 (2008); J. G. Danzl 



5 



322, 231 (2008); J. Deiglmayr et al, Phys. Rev. Lett. 

101, 133004 (2008); F. Lang, K. Winkler, C. Strauss, 

R. Grimm, and J. Hecker Denschlag, Phys. Rev. Lett. 

101, 133005 (2008); S. Ospelkaus et al., Science 327, 

853 (2010); K. K. Ni et al., Nature 464, 1324 (2010). 

J. G. Danzl et al.. Nature Phys. 6, 265 (2010). 
[2] M. H. G. de Miranda et al. Nature Phys. 7, 502 (2011). 
[3] M. A. Baranov, Phys. Rep. 464, 71 (2008); T. Lahaye, 

C. Menotti, L. Santos, M. Lewenstein, and T. Pfau, Rep. 

Prog. Phys. 72, 126401 (2009); L. D. Carr et al. New J. 

Phys. 11, 055049 (2009). 
[4] A.S. Jensen, K. Riisager, D.V. Fedorov and E. Garrido, 

Rev. Mod. Phys. 76 (2004) 215-261. 
[5] F. Ferlaino and R. Grimm, Physics 3, 9 (2010). 
[6] E. Nielsen, D. V. Fedorov, A. S. Jensen, Phys. Rev. A 

56, 3287 (1997). 
[7] F. F. Bellotti et al, J. Phys. B; At. MoL Opt. Phys. 44, 

205302 (2011); F. F. Bellotti et al, Phys. Rev. A 85, 

025601 (2012). 

[8] L Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008); S. Giorgini, L. P. Pitaevskii, and S. 
Stringari, Rev. Mod. Phys. 80, 1215 (2008). 
[9] M. Randeria, J.-M. Duan, and L.-Y. Shieh, Phys. Rev. 
Lett. 62, 981 (1989); Phys. Rev. B 41, 327 (1990). 

[10] A. Pikovski, M. Klawunn, G. V. Shlyapnikov, and L. San- 
tos, Phys. Rev. Lett. 105, 215302 (2010); N. T. Zinner, 
B. Wunsch, D. Pekker, and D.-W. Wang, Phys. Rev. A 
85, 013603 (2012). 

[11] A. C. Potter, E. Berg, D.-W. Wang, B. L Halperin, and 
E. Demler, Phys. Rev. Lett. 105, 220406 (2010). 

[12] S.-M. Shih and D.-W. Wang, Phys. Rev. A 79, 065603 
(2009). 

[13] M. Klawunn, A. Pikovski, and L. Santos, Phys. Rev. A 

82, 044701 (2010). 
[14] A. G. Volosniev, D. V. Fedorov, A. S. Jensen, and N. T. 

Zinner, Phys. Rev. Lett. 106 (2011), 250401. 
[15] N. T. Zinner, J. R. Armstrong, A. G. Volosniev, D. V. 

Fedorov, and A. S. Jensen, ,arXiv:1105. 6264/ 1. 



[16] J. R. Armstrong, N. T. Zinner, D. V. Fedorov, and A. S. 
Jensen, arXiv:1106.2102i^2. 

[17] J. Deslippe, C. D. Spataru, D. Prendergast, and S. G. 
Louie, Nano Lett. 7, 1626 (2007). 

[18] R. Matsunaga, K. Matsuda, and Y. Kanemitsu, Phys. 
Rev. Lett. 106, 037404 (2011). 

[19] D. G. Lidzey et al. Nature 395, 53 (1998). 

[20] X.-Y. Zhu, Q. Yang, and M. Muntwiler, Acc. Chem. Res. 
42, 1779 (2009). 

[21] A. Micheh, G. Pupillo, H. P. Biichler, and P. Zoller, Phys. 
Rev. A 76, 043604 (2007). 

[22] D. DeMille, Phys. Rev. Lett. 88, 067901 (2002). 

[23] A. G. Volosniev et al, J. Phys. B 44, 125301 (2011). 

[24] J. C. Cremon, G. M. Bruun, and S. M. Reimann, Phys. 
Rev. Lett. 105, 255301 (2010). 

[25] J. R. Armstrong, N. T. Zinner, D. V. Fedorov, A. S. 
Jensen, Europhys. Lett. 91, 16001 (2010). 

[26] K. Sun, C. Wu, and S. Das Sarma, Phys. Rev. B 82, 
075105 (2010); Y. Yamaguchi, T. Sogo, T. Ito, and T. 
Miyakawa, Phys. Rev. A 82, 013643 (2010); N. T. Zin- 
ner and G. M. Bruun, Eur. Phys. J. D 65, 133 (2011); 
M. Babadi and E. Demler, Phys. Rev. B 84, 235124 
(2011); L. M. Sieberer and M. A. Baranov, Phys. Rev. A 
84, 063633 (2011); M. M. Parish and F. M. Marchetti, 
arXiv:l 109.2464 i^l; J. K. Block, N. T. Zinner, and G. M. 
Bruun, in preparation. 

[27] C. Mora, O. ParcoUet, and X. Waintal, Phys. Rev. B 76, 
064511 (2007); H. P. Biichler et al, Phys. Rev. Lett. 98, 

060404 (2007); G. E. Astrakharchik, J. Boronat, I. L. 
Kurbakov, and Yu. E. Lozovik, Phys. Rev. Lett. 98, 

060405 (2007). 

[28] Y. Shin, C. H. Schunck, A. Schirotzek, and W. Ketterle, 

Phys. Rev. Lett. 99, 090403 (2007). 
[29] N. Strohmaier et al, Phys. Rev. Lett. 104, 080401 

(2010). 

[30] B. Wunsch et al, Phys. Rev. Lett. 107, 073201 (2011). 
[31] N. T. Zinner et al, Phys. Rev. A 84, 063606 (2011). 



